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We present a general formalism describing a propagation of an arbitrary multiparticle wave packet 
in a one-dimensional multimode chiral channel coupled to an ensemble of emitters which are dis- 
tributed at arbitrary positions. The formalism is based on a direct and exact resummation of dia- 
grammatic series for the multiparticle scattering matrix. It is complimentary to the Bethe Ansatz 
and to approaches based on equations of motion, and it reveals a simple and transparent structure 
of scattering states. In particular, we demonstrate how this formalism works on various examples, 
including scattering of one- and two-photon states off two- and three-level emitters, off an array of 
emitters as well as scattering of coherent light. We argue that this formalism can be constructively 

^Nj used for study of scattering of an arbitrary initial photonic state off emitters with arbitrary degree 

T-H of complexity. 
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I. INTRODUCTION 



_- 1 Recent advance in fabricating quasi one-dimensional (ID) nanostructures has stimulated and motivated extensive 

theoretical works on propagating photons in reduced dimensions. One of the goals of this research is to achieve 

I— I considerable photonic nonlinearities for a single photon pass through the structure which would become a fertile 

i-C platform for variety of applications in future technology and quantum information [l] . The only way nonlinear effects 

may show up in the photonic component of scattering is through interaction with ensemble of emitters. This can 

be achieved by inserting emitters into the ID channel or by side coupling to evanescent modes of the channel. The 

scattering probability of an individual quantum particle of wavelength A in the ID channel is controlled by the ratio 

r^ X^ /A where A is the effective scattering area of the channel. Therefore reducing the effective scattering area makes 

Q^ nonlinear effects feasible in the range of wavelengths close to A. Several physical realizations of such kind of scattering 

'— ' have been recently suggested: tapered optical fibers with atomic ensemble coupled to the evanescent field (currently 

_. 1 tapered region can have optical sub wavelength diameters down to 50 nm) [2 ; hollow optical fibers, "stuffed" with 

<^ cold atoms [3j, photonic crystals [4]. Besides conventional atomic ensembles, various artificial emitters and channels 

^_^ can be used as well. Thus, superconducting q-bits coupled to transmission lines have been engineered to demonstrate 

l/~) two- and three- level emitters (see [5] for extensive recent review). Moreover, photons can be replaced by surface 

^" plasmons coupled to quantum dots ^. 

^^ Having in mind these and other physical realizations we consider the model Hamiltonian which describes interaction 

fT^ of a massless multimode quantized bosonic field propagating unidirectionally in ID geometry with an array of emitters. 

^^ The coupling between bosons and emitters has the standard dipole-like form. The emitters are assumed to have a 

Cn multilevel structure (in particular, in this paper we focus on two- and three-level structures, although our formalism 

^^ allows for a straightforward generalization to more complicated cases as well) . We do not impose any specific details on 

>• the spatial distribution of emitters in the ID channel. However, we distinguish between the case of Dicke arrangement 

(when all emitters are placed in a space region of the size smaller than a typical wavelength) and the case of remotely 

distributed emitters. The problem we focus on is of the scattering type: an initial state consists of a direct product 

j^ of the photonic state and the ground state of an atomic ensemble. Typically, this state is an eigenstate of the free 

(noninteracting) Hamiltonian. An interaction between light and matter is assumed to be switched on adiabatically in 

the infinite past and switched off adiabatically in the infinite future. The adiabaticity parameter controls the so-called 

on-shell condition which plays an important role in the scattering matrix approach (more precisely, the adiabaticity 

parameter defines an uncertainty in the energy conservation condition). The outgoing state is obtained from the 

incoming state by an application of a multiparticle scattering matrix. The main goal is to evaluate the scattering 

matrix, and this is the primary goal of our paper. 

We provide a generic solution of the problem (see Eq. (Is]) below) under the assumptions specified above. It can be 
further extended to the case of multiple emitters. The coupling constant, detuning, and the level structure may vary 
from emitter to emitter. Arbitrary initial condition can also be studied, however we specifically focus here on the 
Fock and coherent states. We found that the multiparticle scattering matrix generically contains many-body bound 
states in the photonic sector, while in the atomic sector it has a projector-like structure. The result of the scattering 
consists in preparing the ensemble of scatterers in a specific linear combination of levels, which is nothing else but the 
dark state, i.e. the state which does not emit. 

Previous studies in eighties [7],[H],[S], nineties (TU], and OO's [H] , [12] , [IS] , [H] , [15] , [E] have already clarified a number 



of theoretical questions related to propagating of photons in ID geometry. Many of these studies have been based 
on the property of integrability. Our analysis here extends these studies in various directions: we consider arrays of 
emitters with multilevel structures and arbitrary coupling constants, and treat the problem in full generality using 
the methodology of the scattering formalism. Our approach does not require any apriori knowledge about system's 
integrability. Nevertheless, our results agree with those obtained by integrability methods, whenever the latter are 
available. 

In the next section we first present our general result with outline of its derivation. In the following sections we 
demonstrate its application on specific examples of light scattering off two- and three- level systems. We consider 
separately the scattering of Fock states with a well-defined number of photons N and the scattering of coherent light. 
The latter study serves the two purposes: first, the scattering matrix in the coherent state basis describes a scattering 
problem of coherent light itself, and, second, it can be considered as a generating functional for scattering problems 
in sectors with a well-defined number of photons. 

II. GENERAL SOLUTION 

A. Main result 

We consider the following Hamiltonian 

H= dvva\v)a{v)+^eiP, + Y^ dvgj{o) {v)S- + a{v)S+} = Hq + a^^Sg,, (1) 



j=i' 



where the spectrum of photons is chiral, that is it contains a single branch w = -\-v corresponding to a unidirectional 



propagation of right-moving photons created by a^{v). Diagonal Pi and off-diagonal S^ = (S*- )^ operators act in 



the Hilbert space of the atomic system, and we use the notations a^ = a+, a = a_, a = ±, a = —a. Here M is 
the number of atoms (emitters) which arc localized in a small spatial region close to a; = 0. The precise form of 
operators 5" is determined by the level structure of an emitter. Thus, for the two-level system in the rotating wave 

approximation (RWA) they read Sj = ^, where cr° are the Pauli matrices. In case of three-level emitters the form of 
operators 3°" depends on a specific type of levels' structure, being either of A, V, or S-type, which will be introduced 
below. Moreover, eventual counterrotating terms may be also included in these generic notations (in particular, for 
the two- level system beyond the RWA we would then have 5+ = 5^ = cr^/2). 

The main assumptions we make by writing up the Hamiltonian M are the following: 1) a linear, unidirectional 
nature of the photon's spectrum and an absence of its lower and upper bounds; 2) the interaction is linear in Cq,; 3) 
a coupling constant (absorbed in Sa) is independent of momentum (energy) v. 

Generically, the scattering matrix is defined as a limit t — > — oo, t' -^ -t-oo of the unitary evolution exponent U{t, t') = 

Texp(— i Jj H(T)dT) and contains different scattering channels. The nontrivial part of scattering is contained in the 
so-called T-matrix, which is related to the scattering matrix S via the expression 

S^\~2m5{E,^-Eout)T{E,^), (2) 

Ein and -Eout being the energies of incoming and outgoing states. 

Under the assumptions we stated above we show that the T-matrix in the A^-particle sector equals 

1 (lu) ^ Ltq G(aQ.j,yjOaj jG . . . G(act2ivi^2iv'-'a2iv j^^o -i (3j 

where the summation over the set of a^ = ± and integration over the set of Vi is implicitly assumed. The hats over 
the operators a mean that they can only be contracted with external operators (that is, operators creating incoming 
and outgoing states). The operators a are effectively normal-ordered in (|3|, but it is still necessary to account for 
their commutation relations with a^ a appearing in the bare Gq ^ = a; — Hq and dressed G~^ Green's functions. In 
turn, the latter amounts to 

G-\w) ^Gq^{u)-Y. = uj-Ho + inS+S-. (4) 

We state that if the self-energy S = —inS+S- does not have a zero eigenvalue, then T(a;) identically vanishes when 
being put on shell (w = Ein = Eout)'- since both Ein and Eout are the eigenstates of Hq, then Gq |os = 0, while G|os 
has finite eigenvalues in all atomic states. To better clarify this important statement we note that the combination 
Gq^G is nothing else, but the projector onto the atomic states with zero broadening (dark states). Therefore the first 
task is to evaluate this building block of (Is]). 



We also note that from the diagrammatic point of view, the class of models we study here does not allow for 
diagrams with intersection and overlapping of photonic lines, what follows from the spectrum linearity and the 
causality imposed by the absence of backscattering for propagating chiral modes. It is naturally possible to extend 
a diagrammatic approach beyond this class, but then the non-crossing approximation no longer provides an exact 
solution, and one should expect effects associated with vertex corrections to the vertices V. 

B. Derivation of the main result 

Below we derive our main expression ([3]). First we remind several basic facts about the scattering matrix approach 
illustrating them on a simple example of scattering off the two-level system. After that, we outline the main steps of 
our derivation. 

1. Scattering problem 
The main goal of the scattering theory is to calculate the scattering matrix 



S — Texp 



+00 



i / V(t)dt 



— 00 



(5) 



where the interaction term V(t) is evaluated in the interaction picture. The matrix elements of the scattering matrix 
Sn'n = WlSln) are defined in the states of the non-interacting Hamiltonian, which in our case have the following 
form 

I") = 4i •■•4„|0)fck)' Sn^'^h + a-, (6) 



i=l 



\n') ^ al^...aljO)b\<j'), e„, = ^ p, -t- a' - , (7) 



i=l 



where |0)b is the photon vacuum state, and a ~ zL labels the states of the two-level system, and fl is the level splitting. 
By definition the T-matrix is 

T{uj) = V + VG{io)V, (8) 

where the full Green's function is defined as G{u!) = {uj — H + iri)~^. The parameter rj controls the adiabaticity of 
switching the interaction on and off in the far past and far future, respectively. The matrix elements of the S'-matrix 
([5]) and the T-matrix (Is]) are related to each other by the the following equation, {n'\S\n) = {n'\n) — 27ri(5(e„ — 
e„/)(n'|r(£„)|n). Therefore, it is sufficient to calculate the matrix elements of the on-shcU T-matrix (that is, at 

UJ = e„). 

2. Calculation of the T -matrix 

We expand G in ([8| in a series of the interaction V 

G{uj) = Go(cj) + Go(a;)V^Go(cj) + Go{u:)VGn{uj)VGo{uj) + . . . , (9) 

where 

Go(c^) = v-r = fT^ + a-^ • (10) 

UJ - Ho + iri LU - ^ - Hb + ir] uj + ^-^ - Hb + ir] 

is the bare Green's function. The projectors P± — <J±a^ — |(1 ± Uz) map onto the spin states |±), respectively. We 
note that P±azp = o-±P± = 0, and P±(j± = a±P^, and therefore T = V + VGqV + VG^VG^iV + . . .. As F is hnear 
in bosonic operators, we can omit odd powers of Y in this expansion: they do not conserve the number of photons 
and will vanish in the calculation of matrix elements. Then T = W + WGqW -I- . . ., where W — VGqV. 



Using the properties of P±, one can show that only diagonal elements T±± 

1 



T++iu) 



g^a^[ 



±|T|±) in spin space are nonzero 

(11) 



T__H 




(12) 



Here the integration over the frequencies {vi^vA is implicitly assumed 

n 

2 

the energy of photons, Wf, = ^^ ki = YliiPi- Then, we can rewrite 



Let us put r++ and T partially on-shell introducing u) = ui, + ^ in Jll|, and w = cjb — -j in (12), where w;, is 



T++{^b) 



T^-H) 



g^au[ 



1 



,t 



LOly 



VL — Hb + irj 
1 



i^b 



5^< 



1 



Hh 



IT] 



LOb 



+ 9^ 



ojb — rj — Hb + irj 
1 



^ uJb — ft — Hb 



IT] 



1 



1 



IT] 



■Ub ■ 



^b — Hb 



17] 



1 

n- Hh 



1 



IT] 



■LOb 



n- Hh 



ITj 



As we are interested in evaluation of average values in the bosonic eigenstates (a^^ ...ap-^T±±a\ 



(13) 



(14) 



it is 



natural to exploit the Wick's theorem, implying (a^/ajj,) = 5{y — v') and {ala^') — 0. We also note the following 
intertwining properties of bosonic operators 

1 1 + 1 1 



.t. 



.t 



' w — Hb — V uj — Hh 



u! — Hb + v' UJ — Hb 



(15) 



An application of the Wick's theorem implies that we have to contract pairwise the operators a^i and aj, . In order 
to do that, it necessary to move a^i to the right, commuting it by means of (15) with all propagators appearing in 
between its initial position and the position of a\,. There are only the two possibilities for such a contraction: 1) 
a^ is contracted with an adjacent operator a\, standing to the right from it; 2) a^, is contracted with some external 
operator a\,. In the diagrammatic representation this means that only diagrams with non-crossing lines of contraction 
are allowed. 

The reason for this restriction is the following. If a^i is reshuffled with more than one propagator, and after that 
it is contracted to some internal operator, we obtain an integral over v' with more than one pole lying in the same 
half-plane. Such an integral identically vanishes. The operator a^' can then be loosely reshuffled to the right end, 
which provides the second possibility. A contraction of the two adjacent operators yields 



LO - Hh 



17J 



-4) = ( / di^diy' 



S{iy 



UJ — Hb — v' 



irj 



(16) 



This generates the self-energy insertion E. 

Let us denote by a^i and a^^ the operators which are contracted to the external operators. In the A^-photon sector 
we have the number N of both species, moreover for the two-level system they alternate (due to u^ = cr^ = 0), that 
is after a^/ we must have a|j, not a^n. All self-energy insertions between such pairs can be resummed, and we obtain 
for r+_|_ and T the following result in the A^-particle sector 



J++ yujb) - 



^2N 



{uJb - Hb 
1 



-iij) 



1 



UJb — Hb 



ivrg^ 



1 



'^ujb + n-Hb 



irj 



1 



UJb 



Hb 



1 



UJb- Hb + iirg'^ 



,t 



UJb 



n- Hb 



irj 



1 



UJb — Hb + ing' 



;{ujb - Hb + irj), 



(17) 



1 1 

V 

Wf, — fi — Hi, + iTrg^ Wb — f^ — -ffb + iirg^ 

I , \ 1 



In these relations operators a and a^ effectively commute wit h each other, and one has only to take carefully into 



account permutations of a{a^) with the propagators using (15). Note that T^^{ijJii) vanishes, when it is put on-shell, 
because of the factors (a;^ — Hi, + ii]) in the beginning and in the end of the corresponding expression. 

Let us now consider T_j{uji,). Moving all propagators to the left and putting them on-shell, we obtain 



T 



(N) _ 2N 



1 



^1 - a All- + i-q All- + M2 - a Yh=i ^j^' + A*w-i - " 

1 1 . . 

a], ...a], a,,' ■ ■ ■ a,,' , (19) 

Ml UN P-M Ml ' y ^ 



where Aai = iii — fit' and a = fl — iirg^ . 

A direct generalization of this derivation to cases with a more complicated level structure leads to our main equation 

III. EXAMPLES 

Here we consider several important examples and compute the scattering matrix for the cases of a single emitter 
coupled to a ID waveguide within and beyond the RWA. We focus on the two-level and three-level cases. 

A. One emitter in the RWA 

For pedagogical reasons, we repeat here the derivation of ( 19 ) from the general expression p). 

In the case of a single atom interacting with a ID field in the RWA we define e± = dz^, P± = a±Uzs^, and 

S± = ga± = g^4^- We find then S = -inS+S^ = -iTTg^P+ and G(w) = P+{uj -^-Hb + ing^)-^ -I- F_ (a; -I- f - 
Hb + if])^^ = P+G+ + P-G-. Therefore, the main building blocks of the equation for the T-matrix are given by 

Go 'G = {u~^-Hb + zri) ,, ^+ . ^ + P-, (20) 

^ oj — ^ — Hi, + i-Kg'^ 

GG,' = ^±^—^(u:-^-Hb + ^v)+P-. (21) 

Putting them on-shell we find 

Lro G|os = GGq |os — P-, (22) 

and therefore 

T{uj) = 527Vp_(a+,,,a_)P+G+(a_,,,a+)P_G_...P_G_(a+.,,„_,f7_)P+G+(a_,,,„a+)P_ 

= P_52^at^G+<G_...G_at^„_^G+<„. (23) 

Introducing w = W(, — ^, we can rewrite G± as G+ = (cj^ — Hb — ft + mg"^)^^ and G_ = (w — Hb + iij)^^ ■ 

B. M emitters in the RWA (Dicke case) 

In the case of many emitters confined in a small (compared to the typical wavelength) region of space we have 
Em = mfl, Pfn = |?7i)(m|, while S± — > gS± . The total effective spin is conserved, SaSa = l{l + 1). If the system is 



initially in the ground state, that is 5*2 = — M/2, then it is sufficient to consider the representation with the largest 
weight I = M/2, to which the ground state belongs. In this case we find 

I 
E = -tirg^S+S. = -ing'{Sl - Sl + S,) = -z^ ^ PmW + 1) - m(m - 1)], (24) 

?n— — ^ 

and therefore for the dressed Green's function we find 

' P 

^ ' ~ 2^ Lj - Hh-Tnn + iTrg'^{l + m){l-m + l)' ^ ' 

The main building block is therefore equal to 

G^^GU = GG^'los = P-i = \- l)H\. (26) 

In the one-photon sector we immediately get 

T^^Huj) - g^P^i{alS-)P^i+,G-i+,{a,,S+)P-i 

Using the action of the collective spin operators on the states of emitters 

S±\m) = .yH T m){l ±m + l)\m±l), (28) 

we obtain that |(— /|5_| — / + l)p = 21. Introducing oj — uji, — fll = uji, — VtM/2 as well as um — ^ ~ iirg'^AI, and 
omitting P_;, we obtain T^^'[uji,) — g'^Mal {ujf, — Hi, — aM)~^o-ii2- Putting this expression on-shell results in 

T^'^^g'M ^ ala,,. (29) 

In the two-photon sector 

T'-^\uj) = g^P^i{al^S-)P-l + iG-l + i{aa2,^,^Sa2)Pm2Grn2{aa3„3Sa3)P-l + lG-l + i{a,y^S+)P-l 

= g^P-i\{-l\S-\ - I + l)\^-l + l|5g, \m2) (m2|5a3 | - / + 1) 

t 1 1 

^"^ LU - Hb + {I - 1)0 + i2TTg^l^°'^''^ Lo-Hb- m2n + ing^H + m2){l - rii2 + 1) 

"" '''''''' Lo - Hb + {I -m + i2^gH''''' ■ ^^^^ 



Using ( 28 1 we obtain 



and 



|(-/|^_| - I + 1)|2(-/ + l|5aJm2)(TO2|^a3l - ^ + 1) 
= {'^I){5a2-^ni2-l^^ + 5a2 + ^m2-l+2\/'2'{2l - I)){5a3 + 5„i2.-l^l + ^a3~^m2-l+2\/'2'{2l - 1)) 
= 2l[2l5a2-6a3 + ^m2-l + "^{^l ~ '^)^a2 + ^a^-^m2-l+'2\ (31) 



r(2)(a;) - ^-'ff'4/V,^^_^^^^^_^^^^^2^^2;«-2^_^^ + ;^ + ,^ 



1 
^ ''''^ a; - iJb + (/ - 1)17 + i2T:gH°"'^ 






P.igHl{2l l)at ,. ^ , „ \,^ , ,,_„„ < 1 



w - i/b + (/ - 1)0 + i27r.g2; "^ uj - H^ + [l - 2)Q + iTTg^2{2l - 1) 
"" ''''^c^-iffc + (/-l)0 + z27r52?''''^- (^2) 



Introducing aM-i = f^ — mg'^{M — 1) and omitting P-i, we arrive at 



T(2)(..,) = 54^/24^ 



1 



1 



.t 



1 



Wh- Hh~ um iOb- Hb + irj ^ ojh - Ht - um 



+ g^M{M - 1)4^ 



1 



.t 



1 



1 



Wh-Hb-aM ^ ^{ujb- Hb) - aM-i ' ^b - Hb ~ aM 



(33) 



Putting T'^2' on-shell, we obtain (note relabeling 1^2 O ly-s in the first term which is allowed because of the symmetric 
nature of the photonic wavefunction) 



1 



1 



r(2) = gH'I- 

Vl — aM Vl +1^2 — K3 - OiM 

In the one-photon sector we have 



M 



M ~1 



Vl - K3 + il-j ^{vi + V2)- OiM- 



^vi^v^^i'a.^i'i- 



+ _A / T(l)t\_x ^^g^ / t t\_A ^^g^ 

vi — aM pi — aM 



and therefore the scattering matrix is 



S, 



'pifei — "piki 



2mt 



pifci — Opifci 



1-2 



iirMg'^ 



Pi - CtM , 
Let us consider the scattering matrix in the two-photon sector. First, we consider 

"Pl+P2,fel+fe2 \^P2'^Pl^ ^ki^k^l ''p2fe2^Plfcl "1" "pifci'-p2fe2 "•" ''p2'£l ''Pi fe2 ' "pi fe2 ^P2fcl • 

Second, 



'5pi+P2,fel+fe2(ap2°^Pl^ Ofcl^fca^ ~ '^Pl+P2,fel-|-fc2g -^^ 3" 



1 



vi - aM Vi - OLM 



M 



- iMnSi,, U-, + 



Af- 1 



and therefore 



where 



V1-V3 2(^1 + '^2)- "A/- 1 

''^\^p\U\^P2V2 ' '^P11^2^P2^1 )\^kil^4^k2l'3 ' ^kil's^k2l'4 ) J 



^PiP2;kik2 — ^pikiJp2k2 ~r i>p2kiJpik2 ^ * 'pip2,fcife2 ' 



rj 



(2) 



PlP2,kik2 



2'K5p^+p^^ki+k29^M 



M Af - 1 

+ T7 ^ 

i^l - i'a 2(^1 + ^2) - "A/^.! 

^\^pil^l0p2U2 ~r '^pil'2"P2l^l)\"kiUi0k2U3 + "fcj jyg Ofcj I/4 j • 



i/1 - aM Vi - aM 



(34) 



(35) 



(36) 



(37) 



(38) 
(39) 



(40) 



Defining E ~ pi + p2 — ki + k2, A — ^{ki — fc2), and A' = i(pi - P2) we find that ki^2 = f ± A, pi^2 = f ± ^'• 
Therefore the T-matrix is 



-(2) 



87r(5p,+p2,fci+fe2g ^'^(l" - "m) 



lp.P2Mk2 [(I _ ^^^^)2 _ A2][(| _ ^^,^)2 _ A,2] 1 ^ ^ i^ 



ttM 



2 ^ Q^Af-1 



87r(5pj+P2,fei+fe2g'*-^(T ^ "Af ) 



2 



[(f - C^mY - A2][(f - aM)2 - A'2] f - CM-l ' 



(41) 



which is in agreement with the result of Ref. |16j . The emergence of the poles in (41) involving more than one 
individual photon's energy refiects a formation of the photonic bound state. 



C. One emitter beyond the RWA 

For the case of an emitter interacting with the field beyond the RWA we have 5*+ = .gcr+ +.g'cr_ and 5_ = gu- +g'a^ . 
Calculating the self-energy we find E = —ittS+S- = —i7:g^P+ — iir^g'YP-, and therefore Gq'^G\os ~ GGq"^|os = 0. 
This implies that the T-matrix identically vanishes. 

Thus, we observe that the non-RWA model is not analytically connected to the RWA one. The nonanalyticity is 
hidden in the noncommutativity of the limits 77 — >■ (which should be taken first for the non-RWA model) and g' — > 0: 
the corresponding eigenvalue 77/(7r(g')^ -I- 77) of the projector onto the emitter's groundstate is sensitive to the order 
of these limits, and a smooth crossover between the non-RWA and the RWA models would be only possible for the 
finite value of 77. 



D. Three-level system, A-scheme 



For the three-level scheme of A-type (see Fig.[l]) we have S+ = g3i|3)(l| -f 532|3)(2|, and S- = g3i|l)(3| -f532|2)(3|. 
The self-energy is S = — i7r(g|]^ + 932)^^ = —i^^yPs, hence the dressed Green's function acquires the following form 



G 



Pi 



P?. 



Ph 



oj - ei - Hf, + ii] LO - €2 - Hb + irj uj - e^ - Hi, + iirg- 



PiGi 



(42) 




I2> 



FIG. 1: Three-level A-scheme 



Therefore the building block is 

Gq G|os — GGq Ios = ^1 + P2- 

The T-matrix in the single-photon sector is 

tW(w) = {P,+P2){alS-)P3G3{a,,S+){Pi+P2) 

= {53'lA+53'2^2+53ig32(|2)(l| + |l)(2|)}at^ 



UJ - 63 - Hb + iirg 



and therefore 



i'pk — ^pk 27ri(73]^dp^: 



Pi 



-2'!Tig3ig32 



Sk 



p + ei - 63 + i7rg2 
|2)(1| 



k+ei,p+€2 ' 



2Trig^2^pk 

Sp+ei,k+i2 



P7 



2 "t'2> 



P + £2 - £3 + «7I"g2 

|1)(2| 
P + ei - £3 + ing'^ 



' P + £2 - £3 + i7rg2 
In the two-photon sector of scattering we obtain (cf. Ref. |13) ) 

T^^\uj) = (Pi + P2)(4i^-)^3G3«^+)(^iG'i + P2G2){alS^)P3G3{a,,S+){Pi + P2) 
= {9iiPi+9l2P2+93i932{\2){l\ + |1)(2|)} alG3a,,i9l,G^ + gl,G2)alG3a,,. 



(43) 



(44) 



(45) 



(46) 



E. Three-level system, \^-scheme 

In the case of three-level V-scheme (see Fig.[2| the S± operators are 5*+ = 53i|3)(l|+52i|2)(l|, and S- = .g3i|l)(3| + 
52i|l)(2|. The corresponding self-energy is therefore E = —iTigliPj, — ing2iP2 — i7r(72i53i{|2)(3| + |3)(2|), while the 
dressed Green's function is G — Pi{uj — ei — Hb + iri)^^ + J2a=2 3 Pai^ — Hb — Aa)~^, where 



P2 = ^2'2^2+e23^3+e2263|2)(3|+6362|3)(2|, 

^3 = ^,33^3 + ^32^2 + e32e33|2)(3| + 6362|3)(2|, 

rh rh 

C22 = C33 



COS-, C23 



-62 = sin-. 



•^2 — ^3 _„„ J, £2 ^ £3 ■ 321 ~ 531 -^2 — A3 



COS = 



A2 + A3 



2 

£2 + £3 



sin0 = -I7r52i53i, 



.921 



9I1 



2 2 2 

(A-A2)(A-A3) = (\-e2+mgl{){\-ej, + iTTgl{)+'!T'^gl^gl^ 



A2A3 = (£2 - J7r.g2i)(e3 - ^-^9:^) + tt 52i331 



(47) 



I2> 




I3> 



FIG. 2: Three-level V-scheme 



Using these relations we observe that our basic building block Gq ^G|os = GGq ^|os = Pi, and therefore 



Pi&l 



1+COS0 2 I l-cos0 „2 I „ „ „• , 
2 321 + 2 531 +g2ig31&in' 

LO - Hh- X2 



+ 



l — costj> 2 
2 921 



i±f^5|i-32iff3isin( 



w - iJf, - A3 

Combining the terms in square brackets, we obtain for the single-particle scattering matrix 

i'^92i{P + ei - £3) + i'^9li{P + ei - £2) 



^pk 



^pk 



(p + ei - £2 + ™9ii){P + ei - £3 + iTT^^i) + T^'^9ii9ii 
J (p + £1 - £2 - iT^92i)ip + ei ~ £3 - ^TTffli) -f TT^gJigli 
*''' (p + £1 - £2 + iT^9li){p + £1 - £3 + iTrglJ + ii'^gligli ' 



(48) 



(49) 



F. Three-level system, E-scheme 



In the case of S-scheme (see Fig. |3| the role of S± operators is played by S+ — g32|3)(2| + g2i|2)(l|, and S- 
532|2)(3| -|- g2i|l)(2|. The corresponding self-energy is S = —i'^932P3 ~ '^'^92iP'2 and therefore G — Pi{uj ~ ei — Hb 
iff)-^ + P2iuj - £2 - Hb + iT^gli)-^ + PsiuJ -e3^Hb + ingl^)^^ = P^Gi. 

The building block is Gq^G\os = GG q^\os = Pi, and therefore 



T(i)(a;) = Pi(4^5_)P2G2(a.,5+)Pi = P^gl.a 



w - £2 + iT^gii 



(50) 



10 



g 



23 



g 



12 



|3> 
2> 

1> 



FIG. 3: Three-level E-scheme 



while 



Spk — 5 



pk 



1 - 2TTi- 



ali 



P + ei-f-2+ 1-^921 



^S. 



pk 



£2 - 17^921 



p + ei - £2 + tTr92i 



(51) 



We note that for E-scheme the one-photon scattering is the same as in the case of the two-level model, the third level 
being inefficient. 

The results for the one-photon scattering of subsections D, E, and F agree with the ones derived in Ref. [12 . 

IV. A MODEL WITH SEVERAL EMITTERS 

Here we apply the same formalism to the case of several emitters. We will distinguish the case of distributed 
emitters from the case of concentrated system (located at the same point) . Moreover we will mostly focus on the case 
of two-level systems and discuss the cases of single- and two-photon scattering separately. We show that our main 
result works even when different emitters have different coupling constants. These calculations explicitly prove one of 
our main statements about an independent character of scattering of unidirectional photons in distributed systems. 



The case of two atoms 



di^i^a\iy)a{iy) + ^a^^'' + ^af'^ + j dv{aHv)S-{y) + a{y)S+{v)} 



To be specific we will focus on the following model Hamiltonian with two emitters 
H = 

where the following combinations play now the role of spin operators 

S+{v) = 9i^t^e^'''''+92afe^''^\ 



(52) 

(53) 
(54) 



It is important to remark that this form of coupling is effectively energy (momentum) dependent, and for ri ^ r2 it 
does not fulfill one of the assumptions we used to derive (Is]). Therefore, an extension of (Is]) is required, which will be 
done below. 

It is convenient to introduce the following projection operators Pab ~ Pa Pjj \ a,b = ±, such that ^^ b-^ab = 1- 
We find then 



^^(l) + ^^(2) = ^(p(l)_p(l))(p(2)+p(2))+^(p(l)+p(l))(p(2)_p(2)) 



Qi + Q2 



P+ 



Ql -^2 



ip+--p-+) 



Ql +^2 



p_ 



(55) 
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It follows that the bare Green's function is 

Go(w) 



R 



+H 



+ 



uj-Hb- ^-^^ +iv io-Hb- ^^J^ + If] 

^± + ^ 

uj-Hb+ ^-^^ +111 Lo-Hb + ^^i±^ + If] 



Pab^ab 



(56) 



B. Calculation of the self-energy and scattering matrices 



There is the only allowed diagram for the self-energy 

I](^)(w) = ai,Sj^yPabGab{!^)alS^^ = S+„PabGab{oJ - l^)S^y 

Pab^ab 



gi<7^+'>e"''-'+92<J^+'^e"''-^ 



giCTi^^e-*'^''i+.g2fTi^^e"*'"'^ 



giai^^e-''"'i+52CTi^^e~"'''' 



= J d,.[glP+bG^b + 9i92e''''-'''-'''^a^-^cT^+^G^^ + g,g,e'''^^^--^^a^^^a^^^G^^ + glPa,+Ga^^] 



(gf + gi)P++ + 9lP+- + 9iP-+ + gig2l{r2.ri)a'ya'f + 5i52/(ri, r2)aV'V^^ 



(1)^(2) 



Jl)^(2) 



(57) 



where f{ri,r2;ujb - Hb) = f{ri,r2) = 2exp[i(wfc - Hb){ri - r2)]8(ri - r2), and uj ^ ujb - ^, ^b being the energy 
of incoming (outgoing) photons. The standard symbol for the Wick contraction is also used. The dressed Green's 
function is therefore 



G(w) = Giujb - Hb) 



P+ 



P_ 



M, 



(58) 



ujb - Hb - rii - V,2 + in{gf + g^) cjb - Hb + irj 

where M. is the matrix to be specified below. 

Now we consider the two cases separately: the case of concentrated system (when positions of two emitters ri and 
r2 coincide) and the case of distributed system when the coordinates of emitters are different. In the former case the 
analysis based on Eq. p| is sufficient, while the latter case requires its extension. 



1. Concentrated case 



For coinciding positions of emitters ri = r2 and we have /12 = /21 = 1 and 

M = I [{G\_)-^ + iTigl] P+_ + [(G"_+)-i + ingl] P_+ + nTg^g2 
1 



(1) (2) , (1) (2) 



detA^ 



[(GO_+)-i + ingl] P+_ + [{G\_)-^ + zngl] P_+ - tTTg,g2 



(1) (2) , (1) (2) 



where 



2//-<0 \-l I „2/^0 \-ll 



detX = {Gl_)-^{G%)-^+z7:[g',iGl_)-^+gi{G%) 

= \ujb-Hb A+ Wfc - i/fc A_ 



A± = - '^(gl +g') ^ l^_^2(g2 ^ g2)2 + (^^ _ ^^)2 _ 2»^(^^ _ ^^)^gl-gl). 



In the limiting case ^\ = ^2 = ^ 

P+- + P- 

2 
P+- - P- 



M = 



1 1 

cjb - _ffb — rj + iTT{g'l + .g|) ojfc - -ffb - fi + i?? 

cos6' 



(1) (2) , (1) (2) 



U!b — Hb - fl + in{gf + p|) ujb - Hb - fl + if] 
1 1 



LOb — Hb — ^ + iTT{gi + .g|) Wb - 7/5 - fi + ir^ 



(59) 



(60) 
(61) 



(62) 
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where tan | = — . Then 

2 91 



^(1) = 



(1) d(2) , n(l) (2) 



G{v,) 



.9i4^^Pi^^+52PW.f 



^Vi'-''V2 '^iyi"'1^2-' 



and the corresponding scattering matrix is 



5pfc — Op]i 



1-2 



iAal + gl) 



p-n + iT:{gl + gl) 



— Spk 



p - n - inigl + gl) 
p-n + in{gf + gf) ' 



gl + gl 



vi-Vt + iTr{gf + .g|) ' 



(63) 



(64) 



2. Distributed system 



For distributed system ri ^ r2 we obtain 



M 



P+- 



+ 



P-, 



u!h — Hb — Oi + ingf LU(, — H(, — 0.2 + JTr^l 

ff \ (1) (2) I ff \ (1) (2) 



*7r5ig2 



{ujb - Hb -fli + iTrgl){uJb - Hb - ^2 + iT^gl) 



= AfiP+_ + M2P-+ + Afi2(T^''(Ti'' + Afaiai'V^''. 



The one-photon T-matrix amounts to 



,1/1 "-"u, 1/2 '-^ + ,1^2-' 



.giCTi^'pl^^e-*"!'^! + g2P^}'^ (t''^^ e-'"^''^ 



G(i/i) La^^^Pl^^e'"^^''! +52Pl^V!^^^e^''^'^^ 



t p dlji'i -^2)+gl{yi -»i) 

and the corresponding scattering matrix reads 

iT^gl{p - ^2) + iT^gliv - ^i) 



Spk — ^pfc ( 1 ~ 2 



(p - ill + i-^gViip -^2 + i-^gl] 



vn 



{p-0^ -mg'i) 



We note that ( 67 1 can be represented as a convolution 



Jpk — 



/ dk'S2-pk'Si-k'k 



(65) 



al,a^2 



(66) 



(67) 



(68) 



of the scattering matrices 5"! and 5*2 on the first and the second emitters, respectively. This property is a consequence 
of the absent backscattering for chiral photons, and it will be later on generalized to the arbitrary iV-photon sector. 
We also make the two following observations: 1) In the absence of backscattering there is no interference between 
counter-propagating waves, and therefore the outgoing state contains no information about positions of emitters, that 
is there is no dependence on ri — r2 in (67). 2) The result (67) would smoothly cross over to (64) on a scale of the 



phononic wavelength. The latter quantity is of the order of an inverse bandwidth, which is effectively set to zero 
in our theory. Therefore, (64) and (67) are not analytically connected with each other in the limit of the vanishing 



distance ri — r2 between the emitters. 



C. General approach: One-photon scattering 



Here we present a general approach aimed at calculation of the scattering matrices for the distributed system. 
Therefore in the following we consider only the model with two emitters located in different positions and having 
different fli and g^. 



H 



diy ly a* {i')a{h') H — 7r<^z 



n 






Vi + V2, 



(69) 
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where Vi — Vi + vj , where i = 1,2. Here Vi = gicr^ J dve^'^^^a{v) = gi(7_l Ai. Let us label the atoms in such a way 
that ri > r2. Our aim is to calculate the T-matrix 

T{u)^V + V Ir^^V, (70) 

where T^ = 14 + V2, in the ground state of the atomic system | 4,4,). In the following we omit irj assuming w — > w + i?y. 



^c.- 


-ffo 


-Vx- 


-V^2 


t 




1 




-'uj- 


-ffo 


-Vx- 


-V^2 


t 




1 




'uj- 


-Ho 


-Vx- 


-V^2 


t 




1 





The first term in (70) can be neglected as it is off-diagonal in spin states. In calculation of the second term 
(Vi + V2) ^_jj ^v +v (^1 '^ ^2) we can retain only (u{ + Uj) in the left factor and {vi + V2) in the right factor, 
respectively, since (t!_| H) = and {H \a_^ ~ 0. Thus we have a sum of four terms 

T = 4- F^^TT TZ^i (71) 

V2 (72) 

V2 (73) 

'^u-H,-V,-V2'''- ^'') 

As a consequence of the RWA we obtain the operators AJ ^ vJ to the left from the resolvent, and the operators 
Ai ~ Vi to the right from the resolvent in each term of this sum. 

Let us now consider term by term expanding them first in Vi, and then in V2. In the following we will use an 
important observation that vi can be only paired with the adjacent u{, while any contraction of V2 and v[ always 
yields zero. More generally, the contraction 

I I. 

v,v] =0, i> j. (75) 

vanishes as it typically implies an integral of the kind 

dv- -— -— = 0. (76) 

{. . . - V + irj) . . . {. . . - V + iri) 

This integral is zero as it has all poles in the upper half-plane, and the exponential function e"'^'''"''^) decays sufficiently 
fast in the lower half-plane for r^ < r j . Therefore we can close the integration contour in the lower half-plane and get 
zero. 



(711 



Let us consider different components (71 1- (74) of the T-matrix in more detail. In the contribution Tn given by 



^''^'''\.-H,-V2^.-H,-V2'''u-H,-V2''''.-H,-V2^--T' ^^^^ 

we can omit the terms in parentheses which are odd in Vi, as they are off-diagonal in the spin states of the i = 1 
atom. Moreover, the operators vi and v\ must alternate. 

If we are exclusively interested in the single-photon scattering, then we should only consider terms with the single 



creation (annihilation) operator in the left (right) side. For this reason we can neglect V2 everywhere in (77 1: an 
expansion in V2 cannot contain V2 in the leftmost position because of the spin state of the second atom, and it cannot 
contain v\ because of the photon state. Note that this argument is implicitly based on the RWA. 
Resumming the remaining series we obtain in the one-photon sector 

where Ei — —iirgi = —iTi. 



Applying similar arguments to T22 (72), we find 



^^ - - L-i/l-E/ - - ^'^^ K-H,-n2-^2 ^^' ^''^ 



-•2 
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where E2 



-iTTgf 



-iFo 



Analogously we find for T12 ( 73 1 that 



Tu 



Vl- 



u) — Ha — V2 w — Hq — V2 



V2- 



(80) 



The terms in parentheses which are even in Vi are omitted; an expansion starts from vi. Once a term containing v\ 
occurs (e.g. ^ wiwjwi), a pairing of adjacent vi and w| should be performed, as there is no any other possibility for 
v\ to be paired in the one-photon sector (the latter implies that v\ can not be either paired to an external vertex). 
This leads to an expression 



'r(i) _ ^t 



Vl- 



LO — Hf) — V2 — ^l UJ — Hq — V2 



-V2. 



(81) 



In the first propagator one can neglect V2 as the corresponding expansion can start neither form V2 nor from v. 
After expanding the second propagator in V2 one keeps only the terms odd in V2, the expansion starting from v. 



Resumming the series containing the power of contraction between V2 and Wj, one transforms (81) into 

1 






1 



1 



-Vl- 



„t 



(jj — Hq — Zjx CJ — Hq to — Hq — Zj2 

It now only remains to pair vi and Wj which results in 



2„2 ^t 



T[:,> = -2mgtg^Al 



1 



OJ — Hq — ill — ^1 



J{ui-Ho)iri-r2) _ 



V2- 



1 



LO — Hq — O2 — S; 



-A. 



(82) 



(83) 



where we have also used that 



uj - ?cr^ LO — ^a^ - ^ 



2"z 



Finally, we consider T21 (74) 



T21 = V. 



OJ-H0-V2 'uj-Hq-V2 



Vl. 



(84) 



(85) 



The terms in parentheses which are even in Vi are omitted; an expansion starts from v[. However, this operator 
cannot be paired to any operator standing to the left from it as follows from (75). Therefore, T21 identically vanishes 



in the one-photon sector, Tjy ~ 0. 

Let us now collect all the terms and calculate the scattering matrix in the one-photon sector 

Ti , r2 _,_ 1 ,„/„__„.^ 1 



Jpk 



Opk 



Opk 



2iS. 



pk 



1 



p — Qi + iTi p — CI2 + *r2 
^^l(p-^!2)+^^2(p-^!l) 



2iTir2e 



-ipr\ 



p - r^i + iVi 



Jp(ri-r2) _ 



p-il2 



{p-ni + iri){p-n2 + ir2) 



Spk 



(p-Qi- ^^l)(p - 1^2 - iT2) 

(p-f7l + ^^l)(p-02 + i^2)■ 



which coincides with the result ( 67 ) from the previous subsection. 



„ipr2 



ir2 



(86) 



D. Two-photon sector 



In order to find an exact expression for a scattering matrix in the two-photon sector it is necessary to classify all 
possible arrangements of external vertices with their eventual renormalization which are allowed by the algebra of 
spin operators. In between vertices one can insert the dressed Green's functions given by Eqs. (58) and (65). 

The diagrams which do not vanish in the two-photon sector are shown in Figs. |4]and[5J The first and the forth 
vertices in these diagrams are not renormalized, which is a consequence of the RWA. This allows us to evaluate the 
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FIG. 4: Diagrams without renormalization of vertices and with renormalization of a single vertex. Sohd Unes correspond to 
the dressed Green's functions, circles denote spin operators, dotted lines are (paired and unpaired) photonic lines. 



FIG. 5: Diagrams with renormalization of double vertices. Notations are the same as before. 



part which is the same for all diagrams 



< (giai^^Pl^^e-^'^i'-i + .g2P!^Vi^)e-"'i''^^)G . . . Ga^^{gia^_^^ P^_^'> e"''"'' + g2P^^^cr^+''e"'^'-') 



aUgiAha'yP'_'>e 



(1) p(2) -ii/iri 



g2M2Pl'a'_'e 



(l)„(2)„-»>.ir2 



g,Ah2P^^'a)_^'e 



WA2)-^u,r^^ 



X ... X a,jMm(r^+^P-^^e"''"' + M2g2Pl'^4'^e''"'' + Mia.gaai'Vl'^e''^^'-^) 



at (giMi^i'^Pi'^e 



(1) p(2) — wiri 



g2M2SiP^''a'_'e 



(l)_(2)„-wir2^ 



X ... X a,,(Mi5i524'^^-''e"'''' + Af2g2Pl''^i'^e'"^''^), 



(87) 



where Si^2 = 1 ~ 27rig^ 2-^^1,27 £^nd Af 1 2 are defined in (k35|. 

Let us first evaluate the diagrams shown in Fig. HI The diagram without vertex renormalization is given by the 
expression 



-* ^jy. ^— .U^ ^^lyn^+.Ur) ^^ly.. ^— .U-^^^Ua ^ + .iyA ^— 



J dj.^O— jj, Lsrdj.O^ i/^(jrCLi,^0^ jj..\srCLj,. O. 



,1^1 ^"1/2 



,1^2 "■ 



1^3 ^ + ,1^3 ^-'"^Z/4 ^ + ,2^4 



P_ 



The first term in ( 88 1 reads 



V J Vl — V-i + IT] 

= {g{Ah{vi)Mi{vi)S2{vi)+g^2Si{i^i)M2{vi)M2{vi) 



vi-vz + irj 



4i«L 01.30^^. 



The second term in ( 88 1 reads 






^b- Hb- ai- Oi2 

X ^ {Nh{u^)S2{vi)e''''^-+'''^'-^ + M2(i.4)e^'^^'''^+'''^"^)<4,a.3«.4 

i!/ — Qfi — a2 

= {M^{vi)Mi{v^)S2{vi)e''^'''~'''^^'''-'-'^ + M2(j/i)5i(i^i)M2(z^4)e'''^^~"^^<''^"'^^^ 
+Af2(z^i)5i(i/i)Mi(i.4)^2(i^4)e*'''^-''^)('^^-'^^)+Afi(i^i)M2(z^4)e-^(''^-'^^)('-^-''^)) 
P— ff^gj t t 
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(89) 



(90) 



(91) 



where ai^2 = ^1,2 ~ *'''5i 2 ^^"^ i? = i^i + 1^2 = '^s + ^4- 

In order to evaluate the diagram with renormalization of the second vertex we find the following vertex correction 
to the second vertex 



(.gifT^a^e''"-^) G [g^a^l^le-^^^^-^) G {g2CT'^^\le-^^^^) 



gtg2a.G^-al^G+^alP\_^'<jl^'e 



t p(l)rr(^^P*''('"l^''2)g^«''2'-i 



Ji^{ri~r2) 



{uJb - Hb - V + iVl){uJb - Hb - f + 1^2 - OLi) 



(92) 



Such vertex correction can only occur in the diagram which has v\ in the first position (from the left). Therefore we 
can effectively replace ujb — Hb — > i^i- Evaluating the integral we obtain 



p|^Vi^^e-*''^''ig2g2(-27ri)e*''i(''i-''=^)Mi(i/2) [l - e^^'^^-^i^^''^^''^) 
P^^'^ Cf^^\'"'^''^ g2e"'^^''^-''^\Si{u2) - 1] 



^ _ pi{i^2-ai){ri-r2) 



(93) 



Now we can evaluate the whole diagram with this vertex correction 



Mi(z/i)[5i(z/2)-l] 



1-e 



i{iy2-ai){ri-r2) 



X (Afi(i/4)5'2(i^4)e'(''^~''^'^('^^~''=^^ +M2(z^4)e'(''^"''*^(''^"''"^) 



P—9h2 

E — ai — a2 



'^l/l'^l/2'-'l'3'^!^4- 



(94) 
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In order to evaluate the diagram with renormalization of the third vertex we find the following vertex correction to 
the third vertex 

gW(ri-r2) 



^(l)p|2)g..3..g^g2^^^ I d„- 



(95) 



{ujb - Hb - ly + 1^3 - a2)iu}b - Hb - v + i-q) 

Such vertex correction can only occur in the diagram which has V2 in the forth position (from the left). Therefore we 
can effectively replace Ub — Hb — > i^a- Evaluating the integral we obtain 



.(1) p(2) ii'3r2 _, '2 



1-e 



= cr^^^pf ^e*''^''^gie*''^(''i"''^) [5*2(^^3) - 1] [l - e*(''3-a2)(n-r2) 
Now we can evaluate the whole diagram with this vertex correction 



i(i'3-a2)(ri-r2) 



(96) 



X [S2{y3)-l] 



_ pi{'^3-a2){ri~r2) 



1-e 



M2{i'i)- 



0,1, Cll„ O-v, CLv 



£/ — ai — a2 
The diagram with renormalization of both the second and the third vertices is given by 



(97) 



J{vi-V3){ri_-r2) 



1-e 



i{v2-ai)(ri-r2) 



2^ _ ^i{iJ3-a2){ri-r2) 
,2„2 



X Mi{v^)[Si{v2) - 1] [^2(^3) - 1] ^2(1.4) 



P-9i9i .t 
E — ai — 012 



'^l/l^I/2^I^3^l^i- 



(98) 



Summarizing the results of Eqs. (91), (94), (97 1, and (98) we obtain the following intermediate expression for T*^^) 



iTr6,,,, [gtMi{iyi)Mi{iy4)S2{i^4) + gtSi{iyi)M2{iyi)M2(jy4)] 
iT^Kusghl [Mi(i/i)M2(i/4) + Mi(i/4)52(;^4)M2(!^i)5'i(i^i)] 

[g4Mi(i.i)Mi(l.4)52(^4) + 52^l('^l)^/2(f^l)M2(i^4)] 



Vl - V3 
9 9 



Vl - 1^3 



[Mi(l.i)M2(z^4) - Ml(i'3)^2(i^3)M2(l^2)^l(i^2)] 6 



-i(i/i-t/3)(ri-r2) 



2 2 
9l92 ^-i(^i-t.3)(ri-r2) 



E — ai — a2 
X {Mi(z.i) + M2{iy2)Si{,y2) + Mi{'y2)[Si{iyi) - 1][1 - e^(-i-"i)('-i-'-^)]} 

X <IVh.[VA} +Ml(j/3)i5(l^3)+M,(l/3)li9(l^4) - llll-e*^''''"""^^''!"' 



{Af2(i^4) + Mi{iy3)S2{iy3) + M2 (1/3) [52(^^4) - 1] [ 



r,-r2)]|^ 



(99) 



which has to be convoluted with al^al^a,y^a^^. The projector P onto the ground state of the atomic system is also 

omitted. 



We see that the dependence on atomic positions is still present in ( 99 1 . Exchanging the dummy frequencies vi 4-> V2 
and t^3 -s-> z/4, when necessary, and exploiting obvious identities 



Nh{y2)[Si{vi) - I] = Mi{v^)[Si{v2) - I]. 

M2{u3)[S2{yA) - I] = M2{Vi)[S2{v3) - l], 

we can cast the position-dependent part of (|99|) to 



(100) 
(101) 



9 9 

9i9i 

Vl -V3 



Ah{vi)M2{Vi) [I - SI{V2)S2{V3)] 



^-i(j^i-i'3)(ri-r2) 



glglMi{vi)[Si{v2) ~ l]A'h{v3)M2{vA)S2{v3)e'^'''-"''>'^'''-'''^ 
glglM2{yi)[S2{v3) - l]Mi(j.i)M2(i.2)5i(i^2)e*^''^-"^)('-^-''^) 

-Mi{v{)M2{v4)[Si{v2) - 1][S2{V3) - iy(E-^^-"2)(r,-r2) ^ 



2 9 

9i9i 



E — ai — a2 



(102) 
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A cancellation of ( 102 ) is achieved by an account of diagrams containing double vertices shown in Fig. [s] The first 
possibility to generate a double vertex is provided by one photon line covering the second and third vertices (upper 
panel of Fig. Isl) . The corresponding double vertex reads 

(1) (2) t / 2 -^(..-.3)n [ , e-('-i--2) 



= -gi92cr+ crL'al^a^^ < g^ 



+5|e"^(''^"'^^)''^ / diy 



{v — vi — ir]){i' — E + ai){i' — t^4 — irj) 



piv{ri~r2) 



O • (1) (2) t 

2mgig2(TX' (J'_' al^ay^ 



X <" f,?p-»('^2-I'3)ri 



5ie 



+526"''"""''"^''" 



{v - vi- iri){v - E + a2){iy - va- irj) 



(103) 



Vl - Ui 



Mi(:^2)Mi(^3)e^(^-"i)(''^-''^) 
Af2(^2)M2(z^3)e*^^""^^^''^"'^^^ 



The whole diagram amounts to 



f Afi(i/i)e"''(''i"''^)(''i~''2) Mi(z^4) 

X [5i(i.2) - l]Af2(^4)^^^^^ + [Si{i^i) - l]M2{v^)^^^ 

+A'h{v^)M2{yA)[Si{y2) - l]Mi(;/3)e*''^^""^^^"^~'^^^ 
+Afi(^2)[^2(^4) - l]'^"^^^ + Nh{v,)[S2{y^) - l]^^^4je 



Vl - Z^3 



Ul - U'i 



+M,{v^)M2{u4)M2{v2)[S2{i^:i) - l]e»('^^-"^)('-i-'-^)} . 



(104) 



The second possibility to generate a double vertex is provided by two photon lines covering the second and the 
third vertices, respectively, and at the same time intersecting each other (lower panel of of Fig. Isl) . The corresponding 
expression reads 



G{gla'•l\^e'^^-' 



(gia'^+\.e'^^')G(gia^l^al^e~''' 
X G (g2ai'^4e-^'''-^) G (g2ofa.,e^''-^^) G (920^:^^ ale'^^^^) 
= 5?52^i'Vi'^a,G__al:,G+_a^G__atG_+a,3G__a|,e*"(''i-"^'e*'^(''i-''^)e-'"^'-i+*"^'^^ 



3„3„-ii.2ri+ii>3r2^(l)^(2) t 

'-'+ "- "1/2 "'^s 



-5i52e 
X I dv I d^i 



r,ii>{ri-r2) pilJ.(rx-r2) 



[v - Vl— iT]){u — E + ai){i' + ^ - E — ir]){^ - E + a2)ip - ^4 - irf) 
Let us first evaluate the integral over /.t 

„i/i(ri-r2) 



(105) 



d/j, 



2m 



{v + ^ — E — iT]){fj, — E + a2){p — V4 — irf) 

^i(E-v)(ri-r2) gi(-E-Q!2)('ri -r2) gii^4(»'i -r2) 



{-V + V2){v'i - v) {v ~ a2){v3 - 0.2) [v " V'i){-V'i + a2) 

\ — p}(v-V3)(ri-r2) 2 — gi{^-'^2)(ri-r2) 



27rze^(^"'')(''i-''^'M2(z^3) 



(106) 



1^ — 1^3 V — a2 

Now we have to perforin integration over u. Note that there is no pole aX v = v^ and v — a2. Collecting the 
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contributions from the poles at v = i^i + ii] and v = E — ai, we cast ( 105 1 to 



X <^ \ 

The whole diagram with this vertex amounts to 

{2mfg\glP^^al^al^a,,a,,Kh{vi)M2{v3)Mi{v2)M2{vA) 



(107) 



vi - 1/3 



E — oil — a2 
alglP — ati«t2"''3a^4 \ -Mi{v2)M2{vi)[Si{vi) - l][52(i'3) - 1 

2 _ pi{E-ai-a2)(ri-r2) 

E — ai — a2 

^[^1(^2) - 1][52(^3) - l]h'h{y{)M2{vA)M2{v2)[l - e»(-2-"2)('-i-'-2)] 
-Nh(v^)M2{vi)[Si{v2) - l][S2{'y3) - l]Mi(i.3)[l - e^(-3-"i)(n-r2)]| _ 



I _ p-i{i^l-i^3)(ri-r2) 



Vi - V3 



+M^{yi)[S2{v3) - l][^i(t'2) - l]M2{vi 



(108) 



Collecting all contributions to T^"^', namely (102 1, (104), (108), and the terms '^ gt^gt from (99), we obtain 



27riT(2) = -- 



27ri 



[glNh{ui)Mi{yA)S2{ui) + gtSi{ui)M2{ui)M2{vi)\ 



2^^glgl { [S,{v,) - l]M2{i^i)^'^^^ + M,{v2)[S2(v,) - if^^""'^ 



Vl - V3 

-M^{y2)M2{yi)[Si{vi) - l][S2{i^3) - 1] 

+Ml{vi)[S2{Vi) - l][Si{v2) - l]M2{Vi) 



Vl - I's 



Vl - V3 

1 



E ~ ai — a2 

-[Siiu2) - l][S2iv3) - l]Miivi)M2{vi)M2{l^2) 

-Mi{vi)M2{vi)[Siiv2) - l][S2iiy3) - l]Mi{v3)} . (109) 

One can observe the dependence on atomic coordinates has disappeared. Transforming T^'^' further, we find 

E -2a2 



27riT(2) 



E-2ai 
Ani 



Slivi)Si{v2)M2{vi)M2{v2)[S2{vi) - l][52(z^3) - 1] 

Mi{vi)Mi{v3)[Si{v2) - l]S2{v^)[Siivi) - l]52(z^3) 



- ;^[^l(^l) - mS2{>^3) - 1][SI{V2) ~ 1][S2{V,) - I]- ^ . 

This expression should be multiplied by al, al a^^a^^. Then the irreducible part of S^"^' equals 



S. 



{2)irred 
PlP2-kik2 



-2TriSp^+p^,k,+k2{ap2apiT'^'^^ 0.1^^2)- 



(110) 



(111) 



Note that if one of the atoms is decoupled (say, 172 =0), then iS'2 = 1, and one recovers from ( 110 1 the already known 
result for the two-photon scattering matrix on a single atom. 

It now remains to evalu ate a reducible contribution to the two-photon scattering. To this end we use an expression 
for T^^' and the (5-part of (102) and find that the corresponding expression factorizes into products of the one-photon 
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scattering matrices 

n{2)red _ 



PiP2,fcifc2 - {SpikiSp2k2 + Sp^k2Sp2k^) {l - 27ri [glMi{p2) + glM2{p2) - 2TTiglglMi{p2)M2{p2)] 
-2m [glMiipi) + glM2{pi) - 2mglglNh{pi)M2{pi)\ 

^{2nif [gtMi{pi)Mi{p2)S2{p2) + gtSi{pi)M2{pi)M2{p2)\ 



{2nifglgl [Ah(pi)M2{p2) + Mi(p2)^2(p2)Af2(pi)^i(pi)] 
{2Trif [gtMi{p2)AUpi)S2ipi) + glSiip2)M2{p2)M2{pi)] 
{2nifglgl [Miip2)M2{pi) + M^{p^)S2{pi)M2{p2)Si{p2)] 



n\^Plki('p2k2 + ''pifc2"p2fcl j i^~'~ 

+ 2(5i(p2) - 1) + 2(52(P2) - 1) + 2[Si{p2) - 1][S2{P2) - 1] 

+2(5i(pi) - 1) + 2{S2{pi) - 1) + 2[Si{pi) - l][S2{pi) - 1] 

+ [^l(pi) - l][Si{p2) - 1]S2{P2) + Si{pi)[S2{pi) ~ 1][S2{P2) - 1] 

+ [S^{pi) - 1][S2{P2) - 1] + [Sl{p2) - 1]S2{P2)[S2{P1) - l]Sl{pi) 

^Si{p2) - l][Sl{pi) - l]S2{pi) + Si{p2)[S2{p2) - l][52(pi) - 1] 

+ [5i(p2) - 1][^2(P1) - 1] + [Slipi) - 1]S2{P1)[S2{P2) - l]Sl{p2)} 

iSp,k,Sp,k2 + Sp,k2Sp2kJSi{pi)S2{pi)S2ip2)Slip2) = SpXS^^X + ^Plk2^P2kl- 



(112) 



Let us now prove that the two-photon scattering matrix on two atoms can be represented as a convolution of the 
two-photon scattering matrices on individual atoms, that is 



S. 



(2) 
PlP2,kik2 



1 



U' ^l.' c(2) 



/ rffcWfc25'2;pip2,fcJ^'5'i.fc'^fe^^fejfe2 



2! / dkj^dk^ \S2.p^k{^2-p2k'^ + ^2;p2k[^2;pik'^ + *'^;pip2.fe'i^ 

,(M c(i) , c(i) c(i) ,,r(2) "i 

^ \'^l;k[ki'^l;k'2k2 ^ '^l;k[k2'-^l;k:^ki ^ '' ' l-k[k'^Mik2 ) 



= {SpikiSp2k2 +Sp,k2Sp2ki)S2{pi)S2{p2)Si{pi)Si{p2) 
+ S2{pi)S2{p2)iTl-p^p2Mk2 + Sl{ki)Si{k2)iT2;p-,p2,kik2 



2 / '^^i'^^'^'2:ptp2,k'ik'J'l-k'^k'^Mi 



(113) 



third terms together yield Sp^I^^^^j^^. We get 



The first term is obviously equal to the reducible contribution (112 1. Let us now show that the second and the 

get 

[2irzglf 



^^Pl+P2-M+k2 

(27rz.gi)2 



47ri 



-{E - 2a^)M^ip,)M,{p2)AMk,)M,ik2)S2iPi)S2iP2) 



Ani 



[E - 2a2)M2ipi)M2{p2)M2iki)M2ik2)Siiki)Siik2) 



1 



{A7rygtgtSp,+p2M+k2iE ~ 2ai){E - 2a2)M2{pi)M2{p2)Mi{ki)Nh{k2) 



X / dk[M2{k[)M2{E - k[)Mi{k[)Mi{E - k[) 
Evaluating the last integral we find 

dk[M2{k[)M2{E - k[)Mi{k[)Mi{E - k[) = 



Ani 



(E-ai- a2){E - 2ai){E - 2a2) 



(114) 



(115) 
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and therefore 



U 



+ 



l+P2-kl+k2 

(27rz.9i)2 



^ /^' [E - 2ai)Miip,)A'h{p2)M,iki)A'hik2)S2ipi)S2iP2) 
Am 



^ . iE~2a2)M2ipi)M2{p2)M2{h)M2{k2)Si{h)Si{k2) 

(27r».9^)^(27r».9|)2 



27ri 



"Pi+P2,k 



M2(pi)M2(p2)Mi(fci)Mi(fc: 

l+fc2 m 



_£/ — Qfl — q;2 



^ PlP2,felfc2' 



(116) 



which reproduces the result contained in Eqs. (110) and (111) 



Note that the sequence of 5*1 and ^2 in the convolution (1131 corresponds to the order in which the right-moving 
photons encounter the atoms 1 and 2 along t he lin e of propagation. The presence of the irreducible part T^^-* in the 



(2) (2) (2) (2) 

two-photon scattering makes the convolution (113) noncommutative 5*2 * 5^ 7^ S*} * S'2 > in contrast to the case 



of single-photon scattering (|68|), where the scattering is commutative, S'2 * 5*} = S*} * 52 



The result (110)- (112) generalizes the corresponding expressions of Ref. [3] to arbitrary coupling strengths gi 7^ (?2- 
Being equipped with the convolution property, one can write down as well an expression for the two-photon scattering 
on M emitters with different gi's. 

The calculations outlined above allows us to make the following general statement: the two-particle scattering 
matrix on an array of distributed emitters is a momentum space convolution of scattering matrices on individual 
emitters. This property has a deep reason: the absence of backscattering for unidirectionally propagating photons. 
Therefore, it is expected to hold also in arbitrary A^-photon sector of scattering. 



COHERENT LIGHT SCATTERING 



In this section we consider scattering of coherent light off a single two-level emitter. The aim of this calculation is 
twofold: at first, we describe the certain physical situation relevant for a discussion of the resonant fluorescence [17] 
in nanostructures, and, at second, we provide a generating functional for a iV-particle scattering matrix. 

We start from introducing operators of photons propagating inside a waveguide of a finite length L 



1 



Ck 



L/2 



2tT J -L/2 



dxa{x)e 



-ikx 



„t 



1 



L/2 



27r J -L/2 



dxo) {x)e'^'' , 



(117) 



where a{x) and a^{x) are the Fourier transforms of the operators a^ and 



,t 



a{x) 



12-K 



dka^e 



ikx 



a\x) = 



77 t —ikx 

dkale 



(118) 



The construction (117) explicitly containing the finite length L is very useful as it allows for a treatment of different 
limits such as a finite photon number limit (N = |Q;fep=const, L — > 00) and a finite photon density limit (7V/L=const 

at TV, L — )• 00). 



The operators (117) obey the commutation relations 



[cfc,cj;] = [ck,al] = 

^ikx 



ap, Cj, 



I gj(fe~p)L/2 _ ^-i{k-p)L/2 



2m 



k — p 



= SA{k-p), 



a{x),c^ 



n-K 



e(-L/2 <x< L/2). 



(119) 
(120) 



where S/^(k — p) is a delta function broadened by A = 2tt/ L and reaching the peak value 1/A at k — p. An expression 
of the operator Ck in terms of ak is given by 



Ck ^ dpapS^ip- k). 



(121) 



As an initial state we choose a coherent state in the mode k 



\<^k 



=-i"»=r/2 



E 






-|0) 



-|afc|V2 



00 ni^ 

a J. 



E 



'nk\ 



\nk) 



(122) 
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This state is normalized to unity, which is ensured by the finite waveguide's length L. It is a superposition of Fock 
states with different numbers of photons n^, and therefore we introduce an operator O which projects on-shell and 
accounts for the energy conservation in each photons' number sector. An outgoing state after the scattering off a 
two-level system is thus obtained by applying the operator O followed by an action of the T-matrix, as prescribed by 
PI, 



\h) = {\-2mTO)\au)^\au)-2^ie-\^''\''^Y. ^ < E 5^" / \X{dq^dp^ 5 [Y^p^-Y^i 

1.1 .t 1 . (4)^ 



.t 



kuk — Hb — a krik — Hi, + irj '^" knk — Hi, — a " n^ ! 



-|0) 



(123) 



|afe)-27rie-l"'=l'/2 J^ 



9 \ "fe/2 "fc (Jyik-n 



L 



^ {uk -n)! 



Yl dq^dp^ j d I ^(Pi -Qi)] iYi ^^^1i ~ ^) 



^i=l 



.t 



1 



1 



Exchanging the order of summations J2'^ =i Sn=i ~ IZnLi '^'n =n ^'^"l shifting Uk — > nj. + n we find 



a 



1 



^ %![., o^Y-" r„. _ ..^ „ ^UY-" r^■-^-^-u,•.l"p^•■■"P"-l[5„_l+q„-p„-a][q„-p„+^,;]'^P"g„-a 



|0). 



OO OO 



3,) = |a.)-2.ze-l-lV2^ ^ r^y"^"^^^;^.+«M)!:52„ ^ ( J] rfg.dp. ) M E(P^ "'?»))( H^Afe - fc) 

n=l nfc=0 ^ ^ 
1 



' n,\ 



.t 



1 



-,t 



"" "^M9i + EL2(9»-pO-«]Er=2fe-pO + *'7]''^^'''"^"-M'?n ^^" -" --1^- -" ^-i"^"' 



-1 + 9n - Pri - a\ [Qn - Pn + ivl Qn - a 



|0) 



|afc) - 2^z ^ i^ — j aW^^ j ( n dq^dp, J (5 ( ^fe " ft) J ( H '^^(ft - fc) 



(124) 



X a' 



ral ...a! 



91 + Z]j=2(ft ^ -P*) ~ "] Ei=2(ft - -Pj) + "/] ^'''" ^" ^ [qn-1 + qn - Pn - a] [qn - Pn + iv] ^" qn - a 



\^k) 



Let us now consider the state 



hi"'') = dpi...dpn dqi... dqnS I ^(ft - ft) j I J| S^iqt - k) 



[qi + E"=2(ft - Pi) - "] [Er=2(ft ~ Pi) + *^] ' " [fti-1 + fti - Pn " "] [?« - Pn + «??] 9n 



(125) 



4, ■•■aLlafc)- 



Performing the inverse Fourier transformation [cf. (118)] 



^' V2n J 



(126) 



and substituting it into (125), we obtain 



h 



(n)\ 



,^ s„/2 dxi... dxn j dpi... dpn I dqi... dq„S ^(pi - qi) \ KQ ^A(ft - k) 



[qi + Z]r=2(ft ^ -P') ^ "] Er=2(ft ^ Pi) + **?] ' " [ft»-l + 5'n - Pn - a] [^n - Pn + «??] ^n - a 



(127) 



Introducing new variables Pj = X]"=i Pi ^i^*! Qi — SiLi ft' such that Pi — Pi ~ Pi+i and g^ = Qi — Qi+i, we cast 
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(1271 to 



|7l"^> = —^^jdxi...dx^jdPi...dPnjdqi...dq^5{Pi^Q^)i\[5^{q,-k)\ 



APi-P2}xi. 



„iPnXn _ 



[Ql - -P2 - a] [Q2 - P2 + i^] " ' [Qn-l - ^n - a] [Qn - Pn+ ^V] Qn " « 

X 0^(2:1) . ..a^(a;„)|afc) 



(128) 



(2 

X I / dF2 



-— ^ dxi... dx„ dqi... dq„ I Jl (5A(gi - fc) I e"^'''' -^ 



o — iP2ixi—X2) 



[Qi - P2 - a][Q2 - P2 + iv] 



dPn 



a 



Qn-1 - Pn - a] [Qn - P71 + iv] 



a\xi) ...a^{xn)\ak). 



Each integral over dPj does not vanish, if one can close the contour of integration in the upper half plane. It is 
only possible for Axj = Xj — Xj^i > 0. Then we get 



dP, 



[Q,-i-Pj-a][Qj-P,+tr^] 



Q{xj — Xj-i)2TTie 



iQiAxi 



i(qj-l-a)Axj __ -^ 

qj-1 - a 



(129) 



and therefore 



oo ^ „ / n \ 

|/3fc) = |ckfc) + ^ (-27rig^Q;fe)" dxi . . . dx„ / d^i . . . dg„ JJ (5a (gi - fe) J B(a;„ > . . . > xi) 



X e 



^'^^-^+■■■+^1--'^ -Q _ _ „t(^^) a^xn)\ak) 



J=2 



= l"fc; 



oo „ 

y [—2TTig'^ak) / dxi . . . (ix„8(x„ > . . . > xi) 



tQ CC \//* "T %i Q ' GL ) .^ IC ' -J- T 

X ( / d(?„^^,5A(9„-fc) n / dq.e^"'"' ""^ ^5A(gj-fc) I at(a;i)...at(x„)|afe), (130) 

9n Q; J \ ^^^ J qj a 



where at = ak/vL. 

It remains to evaluate integrals 






dl^^-^S^'i^-k) ^^^ 



dq. 



Jqj{xj+L/2) _ -iqj{L/2-Xj) 



= r 



1 



-e(.Tj- + L/2) + 



^-j(fe-Q)(i:j+L/2) 



qj + iO+ q.j - ^0^ 



1 



k — a \k — a k — a 

2 



e(-.Tj - L/2) 









fc — a 

I _ (,i(.k~a){L/2-x,) 



(k^a){L/2-xi) ^ ^^ \ 0(^/2 - X,) - r^e{x, - L/2) 



-i{k-a){xj+L/2) _ i(k-a)(L/2-Xj) 



k — a J 

e(-L/2 < Xj < L/2) 

e{xj < ~L/2) 



fc- 



(131) 



and 



iqjXj i(qj-a)Axj + i r iq^Xj + i 

dq, : : <5A(g, - fc) = -6-^"^"^+^ / dq, 5^{qj - fc) 



13 



~ial\Xjj^.\ 



~ial\Xjj^.\ 



^IKXjj^X 



plKXj-^l 



1 — e*('^^")(-^/^^'^^+i) 



9j 
ei-L/2 < Xj+i < L/2) 



^-~i(k^a)(xj + t+L/2) _ ^i(k~a)(L/2^Xj + t) 



e(xj+i < -L/2). 



(132) 
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Summing up the both terms (131) and (132) yields 



e(xj+i - Xj) I dqj (5a (9i - k) 



qj-a 



= e(-L/2 < Xj < Xj+i < L/2) 
+ Q{xj < -L/2 < Xj+i < L/2) 



k — a 
k — a 



I _ ^iik-a)Axj^i 



— i{k—a){xj+L/2) _ i{k—a)Axj^i 



(133) 



Note that there is no contribution to (1331 when both Xj and Xj+i are smaher than —L/2. 



Thus, we fin ally obtain an expression for the state \/3k) which emerges after scattering of the initially prepared 
coherent state ( 122 1 off the two-level system 



\Pk 



"fc 



E 



2mg'^ak 
k — a 



dxi...dx„e*'=(^i+-+^"^ 



e(L/2 >Xn> ■■■> X2>xi> -L/2) Y\ \l 



J(k-a)Axj-ii 



+e(i/2 > x„ > . . . > a;2 > -L/2 > a;^)e-»('=-")(^i+^/2) fi _ ^^{k-o.){x2+L/2)\ Y[(l- e'('=-")^^.+i 
X a\xi) ...a\xn)\ak), 



(134) 



where Xn+i = L/ 2. 

The state ( 134 ) has a remarkable property: when two coordinates in the integrand approach each other, it vanishes. 
This property lies in the origin of antibunching of photons which is conventionally observed in the density-density 
correlation functions [I7],[5]. An occurrence of the two contributions with xi > —L/2 and xi < —L/2 has been 
previously remarked in [8], both being important for a proper normalization of the state (134). 

The explicit expression for the outgoing state ( 134 1 opens a possibility for a study of correlation functions as well 



as photons' statistics in all parametric regimes, which will be a subject of a subsequent publication 



VI. CONCLUSIONS 



We have developed the scattering approach to problems of propagating bosons in one dimensional geometry. We 
have derived the general form of the T-matrix (non-trivial part of the scattering matrix), Eq. (13), under the following 
assumptions about the spectrum of bosons: the spectrum is linear, chiral, and infinite. Our formalism is complimentary 
to the Bethe ansatz solutions [7] and to the traditional approaches based on the equations of motion. 

We have applied the developed formalism to several specific examples including emitters with two- and three-level 
structures. The emitters can be either distributed across the ID channel or concentrated in a tiny region of space. 
We have shown that the scattering results in projecting the state of the emitter onto the specific - dark - state which 
does not emit. 

We have also shown that the one- and two-particle scattering matrices off two emitters can be represented as a 
convolution of scattering matrices corresponding to individual emitters. Hereby the microscopic properties of different 
emitters can vary (i.e. coupling constants to photons, detunings, level structure). We conjecture that this property 
generically holds for multi-emitter arrays as well as for the multiparticle scattering with N > 2 photons. We are going 
to elaborate more on these issues in future studies fTSl . 

The developed approach can be applied to the arbitrary initial state, which can either conserve the particle number 
or not. In the case of number-conserving initial state, we observe the formation of photonic bound states, which is 
reflected in the emergence of a pole of the 5'-matrix involving in its argument an energy of more than one individual 
photon. In the case of the coherent light scattering, we clearly observe the fermionized behavior typical to the 
To nks- Girardeau gas [19] discussed in the photonic context in [20| . 
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